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ABSTRACT 

We revisit the case of magneto-rotational evolution of neutron stars (NSs) with ac- 
counting for changes of the angle x between spin and magnetic axes. This element of 
the evolution of NSs is very important for age estimates and population modeling, but 
usually it is neglected. In the framework of two models of energy losses we demonstrate 
that unless NSs are born with the inclination angle xo very close (within ^ 1°) to the 
position of maximal losses, pulsars with short initial periods quickly reach regions of 
small energy losses (aligned rotators in the case of magneto-dipole losses and orthogo- 
nal rotators - in the case of current losses) without significant spin-down. This means 
that either most of known NSs should be born with relatively long periods close to 
presently observed values, or the initial inclination angles should be large (close to 
90°) for the magnetodipolc model and very small (close to 0°) for the current losses 
model, or that both models cannot be applied for the whole evolutionary track of a 
typical NS. In particular, magnetar candidates should be born with periods close to 
the observed ones within the scope of both models of energy losses. We discuss how 
these considerations can influence population synthesis models and age estimates for 
different types of isolated NSs. Focusing on the model of current losses we illustrate 
our conclusion with evolutionary track reconstructions on the P - sinx plane. We 
conclude, that most probably both existing models of energy losses - magneto-dipole 
and longitudinal current - require serious modifications, on the other hand, estimates 
obtained under the standard assumption of sinx = 1 = const does not have solid 
theoretical ground and should be taken with care. 
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1 INTRODUCTION 

Radio pulsars have been discovered nearly 40 years ago 
jHewish at aLlligeSl ). and very quickly they were recognized 
as magnetized rotating neutron stars spinning down and dis- 
sipating their rotational energy. Despite of all the progress in 
the theory of magneto-rotational evolution of isolated neu- 
tron stars, several important questions still remain unclear. 
The uncertainty concerning the evolution of the angle x be- 
tween the magnetic axis and the spin axis of an isolated 
neutron star is one of them, although the angle x signifi- 
cantly affects the rate of energy release for all mechanisms 
of radio pulsars emission. At that, among the basic parame- 
ters, which determine astrophysical manifestations of radio 
pulsar, this one is the least popular, and its evolution is 
usually neglected. 
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We discuss evolution and estimate ages of isolated neu- 
tron stars taking into account changes of the angle between 
the magnetic and the spin axes. We consider the magne- 
todipole model and the model of pulsar radio emission based 
on current losses in the magnetosphere of a neutron star, 
focusing on the latter. It should be emphasized here that 
the evolution of an inclination angle is a key parameter for 
both models. Unfortunately, many authors either completely 
disregard the evolution of an inclination angle or at best 
consider its result to be identical to the magnetic field de- 
cay (see, for examp le, 'Regim bau fc de Freitas Pachecci|l200ll : 
I Fau c her- Gi g uerefc Kaspi 2003, and references therein). 

In the following section of this paper we consider the 
magnetodipole and current losses models of energy release 
of radio pulsars. There the basic information about the mod- 
els themselves and their main contradictions with observa- 
tions are listed. The expression for the characteristic ages 
of isolated neutron stars with the account of evolution of 
an inclination angle is presented. As a result, the ages of 
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isolated neutron stars estimated by the obtained formulas 
for different initial parameters (spin periods and inclination 
angles) are listed in Tables 1 and 2. In the third section 
we reconstruct possible evolutionary tracks of XDIN RX 
J0720.4-3125 and AXP 1E2259+586 varying parameters and 
assuming the model based on current losses is valid. In the 
forth sections we give a discussion. In particular, as the ac- 
count for changes of x mainly influence population synthesis 
models and age estimates, we discuss these topics. Finally, 
we present our conclusions. 



2 THE ROLE OF THE EVOLUTION OF 
INCLINATION ANGLE IN THE AGE 
ESTIMATES ON THE ACTIVE RADIO 
PULSAR STAGE 

In this section we estimate duration of active radio pulsar 
stage for isolated neutron stars taking into account the evo- 
lution of the angle between the magnetic and the spin axes. 
We consider two models of energy release mechanism: the 
magnetodipole model and current losses model. Everywhere 
in the paper we neglect the possibility that magnetic field de- 
cays, and that there are any non -electromagn e tic effects like 

(ll97-^ '): lMenou et all 

mm). 



According to ([T]), the slow-down rate of a radio pulsar for 
the magnetodipole model can be defined as follows: 



caya, anu liiul Liieie aie any lion-eiecuio 

those d escribed, for examp l e. bvlMacvl I 
l|200ll ): IWassermanl (|2003l ): iRudermanI ( 

The primary point of our calculations is that the ratio 
f{x)/P, where P - is spin period, and fix) = cosx for the 
magnetodipole model and /(x) = sinx for the current losses 
model, is invariant during the evolution of a neutron star on 
the radio pulsar stage. Therefore, in order to estimate the 
characteristic age of an isolated neutron star in scope of one 
of the models we do not need to know the initial inclination 
angle. The inclination angle and the spin period of a neutron 
star at the present time and its initial spin period are enough 
to define the initial inclination angle for each model of energy 
release. The validity and implications of the invariants are 
discussed throughout this paper. 



2.1 Magnetodipole model 

Energy losses of an isolated radio pulsar due to radiation can 
be associated with two different scenarios. The first one is 
based on magnetic dipole braking, and we discuss it in this 
subsection. In this case, assuming a perfect dipole rotating in 
vacuum, the energy loss es can be est imated by the so-called 
magnetodipole formula (|Pacinilll967l ): 
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Here J ~ 2/5MR^ is a moment of inertia of a neutron star, 
X is the angle between the spin and the magnetic axes, and 
Q = 2it/P is the angular velocity of a neutron star. The 
magnetic field is defined as Bo = 2/i/7?^, i.e. this is the field 
on a pole. The magnetodipole losses for ordinary pulsars 
are approximately lO^'^ - 10"^'* erg s~^. A simple estimate 
of a lifetime of a normal radio pulsar with a spin period 
P = 1 s and period derivative P = 10"""^^ s s~^ gives us 
Teh ~ P/2P « lO'' yrs. One can see that this result is only 
a rough approximation; for more accurate estimate of the 
age of a radio pulsar the evolution of the angle between the 
magnetic and the spin axes should be taken into account. 



P = 0.24 X lO^^^B^aP"^ sin^ x- 



(2) 



In the scope of the magnetodipole model the radiation- 
reaction torque is perpendicular to the magnetodipole mo- 
ment of a neutron star. Therefore, the component of angular 
velocity SI cos x parallel to the magnetodipole moment (and 
so per pendicular to the deceleration momen t) remains con- 
stant (IDavis fc GoldsteinlHoTOl : lMichellll"99lh . i.e. the ratio 



cosx 



(3) 



is invariant during the evolution of a neutron star at the 
active radio pulsar stage (since birth of a neutron star and 
up to the moment when it crosses the death line). Thus, the 
angle between the spi n and the magneti c axes during the 
evolution tends to 0° . iDavis fc GoldsteinI (|l970l ) emphasize 
that this model has the following implication: if radio pulsars 
are born as fast rotators and are decelerated only by dipole 
radiation to spin periods close to one second then the initial 
inclination angles must be very close to 90°. 

Here we outline the problem that most of the authors 
disregard as they do not take into account the evolution of 
the angle between the magnetic axis and the spin axis. It 
is often accepted, especially in population synthesis models, 
that radio pulsars are born as fast rotators (Po ^ 10—30 ms). 
As there is no information about the initial distribution of 
radio pulsars in the inclination angle (.Tauris fc Manchester! 
Il998l ). most of the authors assume the distribution in the 
initial inclination angle to be equiprobable, or just ignore 
this question. Eq. @ defines the evolutionary track of a 
newborn radio pulsar with given Po and xo till it reaches 
the death line, which for the magnetodipole model is gener- 
ally appr oximated by the following simple equation (see, for 
example, iBhattacharva et al.lll993 ) : 



(4) 



^ =0.34 X 10^^ Gs"^ 
P^ 

It should be noted here that the death line approximation 
((4]) is valid for a neutron star with a perfect dipole magnetic 
field and with the radius R = 10*' cm, the moment of in- 
ertia J — 10*^ g cm^, and the inclination angle sinx ~ 1, 
i.e. B is calculated out of P under these assumptions. The 
lack of dependence of the death line on the angle x should 
be considered as an important disadvantage of the usual ap- 
proach in the context of the magnetic dipole braking model. 
In the model of current losses such dependence is taken into 
account (see below). Accordingly, the accuracy of age esti- 
mates is often lower for the magnetodipole model than for 
the model based on current losses. As a result of @, for the 
period Pdi, at which neutron star is leaving the radio pulsar 
stage (the death line crossing), we have: 

Pdi = 1.7B?2's, 



here P12 = P/IO'^ G. 
Thus, if 

P.>^, 
cos Xo 

an isolated neutron star becomes an aligned rotator before 
it can reach the death line and leave the stage of radio emis- 
sion. As one can see, this way of evolution should be typical 



(5) 



(6) 
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for most of the radio pulsars with short initial spin periods, 
except for those with the initial angles between the mag- 
netic and the spin axes very close to 90° . For example, let us 
consider an evolutionary track of the neutron star with the 
initial spin period Pq = 20 ms, initial angle between the spin 
and magnetic axes xo = 45° and magnetic field Bo ~ 10^^ G. 
Eq. (|3]) allows us to estimate the period at which the neutron 
star becomes an aligned rotator: Pai ~ 28 ms. 

As a result, for the magnetodipole model we conclude 
the following: if the invariant (|3]) is valid, and radio pulsars 
are born as fast rotators, and their initial angle distribu- 
tion is equiprobable, then most of radio pulsars should be- 
come aligned rotators long before they pass the death line. 
Thus, we should expect the distribution in observable an- 
gle between the magnetic and the spin axes of active radio 
pulsars to have a peak close to 0°. Otherwise, an essential 
part of radio pulsars should initial ly have long spin periods 
(Seskin. Gurevich fc Istominlll993l '). 

Further, we can estimate a lifetime of an isolated neu- 
tron star as an active radio pulsar: 



Tmd[s] 



dP 



Po 



2.4 X lO-iesSjP-i (1 - cos2 Xo(P2)/(p2)) 



,(7) 



where Po and xo are the initial spin period and the angle 
between the spin and the magnetic axes of a radio pulsar, 
respectively, and Pdi is a period at which a radio emission 
turn-off occurs (a pulsar crosses the death line). If an inclina- 
tion angle x reaches 0° before a pulsar crosses the death line, 
the period of the transition to aligned rotators Pai should 
be used instead of Pdi. 

Substituting Eq. Q into ((7)l and integrating, we obtain: 



r„,d = 2 X 10'°Pi 



Po 
cos xo 



In 



'XO 



p2 p- 



'XO 



s.(8) 



Therefore, assuming the initial spin period Po = 20 ms, 
initial angle between the spin and the magnetic axes xo ~ 
89.6° and magnetic field Po = 10^^ G, from Eq. (HI we 
obtain an estimate of the age of an isolated neutron star at 
the moment it crosses the death line as rmd = 226 Myr. If we 
assume xo = 45° , the age of the radio pulsar which enters the 
stage of aligned rotators can be estimated as Tai = 0.16 Myr. 
Table 1 contains the estimates of duration of the radio pulsar 
stage for different parameters (such as magnetic fields, initial 
spin periods and initial inclination angles of radio pulsars). 

It should be noted here that the age of a neutron star 
estimated within the scope of magnetodipole model by ([8]) 
significantly depends on the death line. For example, if the 
coefficient in ([5| is twice less, then the lifetime of a typical 
Crab-like pulsar at the radio pulsar stage estimated by ^ 
is reduced by at least the factor of 4. At that, the smaller 
inclination angles x at which the pulsar crosses the death 
line gives the larger decrease in the estimated lifetimes with 
the coefficient reduction in ([5}. As the spin-down rate of a 
neutron star decreases near the death line (due to the factor 
sin^ X in El- 0i the longest part of a lifetime at the radio 
pulsar stage is the period when a neutron star stays near 
the death line. Thus, as we do not take into account the 
inclination angle x for the death line in the scope of the 
magnetodipole model Eq.([8]) only roughly approximates the 
age of a neutron star. It is most likely that the more accurate 



expression for the death line can reduce the age estimates 
by a factor of a few. 

2.2 Longitudinal current losses model 

In this subsection we estimate ages of radio pulsars assuming 
that the energy losses of an isolated neutron star are associ- 
ated with longitudinal currents in the pulsar magnetosphere. 
In our estimations we take into account the evolution of an 
angle between the magnetic and the spin axes of a neutron 
star. The energy losses of an isolated neutron star estimated 
by magnetodipole formula are of the same order of mag- 
nitude as the so-called current losses - the energy release 
determined by ponderomotive forces of longitudinal electric 
currents, i.e. the currents fiowing along the o pen magnetic 
field lines and closing on the surface of the star (iBeskin et al.l 
[l99^y 

A simple age estimate of a normal radio pulsar for the 
model based on current losses is the same as for the magne- 
todipole model: Tch ~ P/2P « lO'^ yrs. A more accurate es- 
timate of a radio pulsar lifetime for the current losses model 
requires to consider evolution of an angle x between the spin 
and the magnetic axes of a neutron star. The deceleration of 
neutron star rotation in the scope of this model is caused by 
the surface currents that flow across the magnetic held lines 
(jSeskin et al.l[l99^ ). Therefore, the ponderomotive Ampere 
force arises, and the direction of the appropriate deceleration 
torque is opposite to the magnetodipole moment. Thus, the 
component of the angular velocity of a neutron star sin x 
perpendicular to the magnetodipole moment (and so per- 
pendicular to the deceleration moment) remains constant, 
and the ratio 



smx 



(9) 



is invariant during the stage of normal radio pulsar. As we 
can see, in the scope of current losses model the angle x 
between the magnetic and the spin axes evolves to 90°. A 
time of evolution of an inclination angle is of the same order 
of magnitude as the slow -down rate of a neutron star. 

iBeskin et af] (|l993[ l developed a theory in which en- 
ergy losses are associated with longitudinal currents in the 
neutron star magnetosphere. In this approach there is no 
unique formula for P for all phases of a radio pulsar evo- 
lution, also different assumptions about particle escape can 
be made. For the escape model based on the proposal by 
Rudorman & Sutliorlanc| (|l975[ l (hindered escape) and for 
pulsars not very close to the death line, the following equa- 
tion is derived: 



P = 10~"P 



15DlO/7pl/14, 



J/2 



(10) 



Note, that the exact value of the coefficient in this equation 
is not well defined. In this paper we use the value written 
above. So, in the scope of the current losses model, we can 
write the following expression for the time that an isolated 
neutron star lives as a radio pulsar: 



-Tcl S 



dP 



P„ 10-i^Bi^°/ Wi4 (1 _ sin^ ;^0(P2)/(P2)) 



(11) 



Po and Xo are the initial spin period and the initial angle 
between the spin and the magnetic axes of a radio pulsar, 
respectively, and Ptr is a period at which neutron star is 
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Table 1. Age estimates for the magnetodipole model. Po and xo are initial parameters of a neutron star, 
Pdl and Tmd are the spin period and the estimated age of a neutron star on the death line, Pai and Tai are 

tiio sj)iij period and tlic estimated age of a uentrou stai' at tlie moment of transition to the aligned rotators. 
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Table 2. Age estimates for the model based on longitudinal current losses. Pq and xo are initial parameters 
of a neutron star, Pdi and Td are the spin period and the estimated age of a neutron star on the death line, 
Port and Tort are the spin period and the estimated age of a neutron star at the moment of crossing the 
boundary of the region of orthogonal rotators. 
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leaving the radio pulsar stage. To find Ptr let us consider 
two possible subsequent steps of the evolution of an isolated 
neutron star. At first, as a radio pulsar spins down, it may 
cross the death line and pass to the region of the so-called ex- 
tinct radio pulsars. These neutron stars have already ceased 
to radiate in the radio band, but electrodynamic processes 
in their magnetospheres still play a crucial role. In this case 
the critical period Pdi is defined by the death line of radio 
pulsars: 



B 



8/15 



/ \0.3 

(cosx) S. 



(12) 

Using Eq. © , we can find an equation which shows how the 
period of transition of a neutron star to the region of extinct 
radio pulsars depends on the initial parameters: 



sm Xo-p2, 



(13) 



dl ^0.15 

If the model based on the current losses is valid, the 
period derivative (|10p reduces during the evolution of a neu- 
tron star due to the factor cos'''''^ x- So, the longest part of 
the lifetime of an active radio pulsar is the time it spends 
near the death line (unless a radio pulsar passes to the stage 
of orthogonal rotators before it can reach the death line - 
this way of evolution we discuss below). Thus, the age esti- 
mate (|11|) significantly depends on the accuracy of the ex- 
pression (|10|) near the death line. Unfortunately, as is was 
noted above, expression (|10|l adequately defines the period 
derivative only for the radio pulsars far from the death line. 
Thus, the age estimates for the current losses model are very 
approximate and the more accurate expression for the spin- 
down rate of a neutron star near the death line can alter the 
results obtained in (|lip by a factor of a few. 

As it is noted in the previous paragraph, a NS in- 
stead of crossing the death line can follow a different path. 
An alternative evolutionary track lead s an object d i rectly 
to the stage of orthogo nal rotators (|Beskin et al.l 1 19931 : 
iBeskin fc Nokhrinall2004r ). Here a short explanation should 
be done. If an inclination angle of a radio pulsar is close 
enough to 90°, so that 



\7t/2 ~x\< Ro/R, 



(14) 



then particles with opposite charge signs can fiow out of the 
polar cap surface. Here Ro = R{iiR/c)'^"' is the polar cap 
radius. A group of radio pulsars for which the Goldreich- 
Julian charge density changes the sign on the polar cap sur- 
face is called orthogonal rotators. These pulsars are charac- 
terized by fast {P ~ 0.3 s) or superfast (P < 0.1 s) rotation, 
reduced spin- down rate (P ~ lO"^'' - 10"^'' s s"^), weak 
radio emission (Lrad ~ lO'^^ — lO'^^ erg s"'^), and the ex- 
istence of an interpulse. The energy losses for this class of 
objects are reduced by the factor Q.R/c in comparison with 
the normal radio pulsars and are virtually independent of 
the angle between the magnetic and the spin axes, which 
is anyway very cl ose to 90 degrees. Some ca ndidates to this 
class are listed in (|Beskin fc Nokhrina|[2003 '). They have pe- 
riods from ~ 100 up to ^ 500 msec. Not all radio pulsar 
with small P should be orthogonal rotators. It was noted by 
iLorimer et al] (l2004h that isolated pulsars with small period 



derivatives. 



10" 



and shorter periods about few dozens 



of milliseconds can be "disrupted recycled pulsars", but ob- 
jects with longer spin periods require a different explanation 
which is probably related to emission mechanism. 



The boundary of the stage of orthogonal rotators is 
given by (|14|l : 



cosx = {2-kR/cP) 



1/2 



1.45 X 10"^P"" ' 



(15) 



so it is easy to find the following implicit equation for the 
period of the transition to this stage: 
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The death line of radio pulsars (|12|) and the boundary 
of the orthogonal rotators stage (|15p intersect at the point 

Pint « 0.35B?2'"' s, sinxint ~ (l - 6 X IQ-^'PiV)"'^ in the 
P - sin X diagram. Therefore, if 

(1 



sinxo/Po > 



6 X 10-^PiV')°-5 



0.35502*'^ 



(17) 



a neutron star passes through the death line to the region of 
extinct radio pulsars. In this case, the period of transition is 
determined by the death line p2p . so that Ptr = Pdi. Thus, 
assuming the initial spin period Po = 20 ms, the initial angle 
between the magnetic and the spin axes xo = 1°, and the 
magnetic field B — 10^^ G, from Eq. (|lip we can estimate 
the age of a radio pulsar at the point of the transition to 
the region of extinct radio pulsars as ~ 36 Myr. The spin 
period of a neutron star at the point of the transition is 
Pdi ~ 0.88 s, and the inclination angle is Xdi ~ 50°. For the 
same spin period and standard period derivative P = 10~^^ 
ss~^, the characteristic age estimate Teh = P/2P results in 
14 Myr. As we can see, taking into account the evolution of 
the angle between the spin and the magnetic axes within the 
current losses model significantly increases the time that a 
radio pulsar needs to slow down from the spin period Pi to 
the larger spin period P2 (both Pi and P2 are supposed to 
belong to the region of active radio pulsars). 

Otherwise, if the initial parameters of a pulsar do not 
meet the above condition (|17p . the transition to the stage 
of orthogonal rotators takes place, so that Ptr = Port. Note, 
that the transition to the orthogonal rotator stage should 
be typical for most of the radio pulsars with comparatively 
short initial spin periods Po ~ 10 ms (the ratio sinx/P is 
initially very large for these radio pulsars provided that the 
angle x is not too close to 0°). Thus, for the model based 
on current losses the majority of neutron stars pass to the 
orthogonal rotator stage before the transition to the region 
of extinct radio pulsars. Assuming the same initial period 
Po = 20 ms and xo = 45°, we can estimate the age of a 
neutron star passing from the stage of radio pulsars to the 
orthogonal rotator stage as follows: Tort = 1.44 Myr (as- 
suming B = 10^^ G), Tort = 54000 yrs (B = 10^^ G), and 
Tort = 2000 yrs {B = lO" G). The spin period at which 
a radio pulsar crosses the boundary of the orthogonal rota- 
tor stage for all three magnetic field strengths is the same 
Port « 28.2 ms. As one can see, it requires about 1 Myr for 
an active radio pulsar with comparatively large initial an- 
gle X a^nd standard magnetic field to increase its period by 
only 8 ms. Table 2 contains the age estimates for different 
parameters (such as magnetic fields, initial spin periods and 
initial inclination angles of radio pulsars). 

Thus, we show that if the evolution of an inclination an- 
gle X is taken into account, the slow-down rate of an active 
radio pulsar turns out to be several times lower, and a life- 
time of a star at that stage turns out to be much longer than 
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it is commonly assumed. Moreover, as the angle x between 
the spin and the magnetic axes of a pulsar increases, the 
slow-down rate of the star reduces by the factor of (cosx)"", 
where the power m is ~ 1.5 - 2. Fu rthermore, as it is dis- 
cussed in l|Beskin fc Nokhrinalliooi) . a period derivative P 
of neutron stars at the stage of orthogonal rotators becomes 
approximately 10* smaller than for the normal radio pulsars, 
so we should see some weakly radiating old radio pulsars 
with comparatively small spin periods and with the angles 
between the magnetic and the spin axes close enough to 90° . 

After all, if radio pulsars are born as fast rotators and 
their initial angle distribution is equiprobable, for the model 
of current losses the most of radio pulsars should become 
orthogonal rotators long before they can pass through the 
death line. Thus, as opposed to the magnetodipole model, 
we can expect the distribution in present day angles between 
the magnetic and the spin axes of active radio pulsars to have 
its maximum close to 90°. If the angle distribution has no 
such peak, we can assume that either the initial inclination 
angles are very small, or an essential part of observable radio 
pulsars was born with long spin periods. The latter conclu- 
sion is quite probable, and has a long history of d i scussi on in 
the literature (see, for example, iHarding fc Lail l|2006l ) and 
references therein). 



3 RECONSTRUCTION OF EVOLUTIONARY 
TRACKS: CURRENT LOSSES MODEL 

As it is shown above, the evolution of the inclination angle 
of an isolated neutron star significantly affects the estimate 
of its characteristic age through the influence of this param- 
eter on the spin-down rate and the shape of evolutionary 
track of a neutron star. Usually, evolutionary tracks of NSs 
are shown on the P - P plot. However, in the case of evolv- 
ing X it is convenient to use the P - sinx diagram. In this 
section we describe the procedure of reconstruction of evo- 
lutionary tracks for the isolated neutron stars assuming that 
the model based on current losses is valid. Two models of 
particle acceleration region should be considered: the model, 
in which the longitudinal electric field near the neutron-star 
surfa ce hinder the particle escape (|Ruderman fc Sutherland! 
1 19751 ) and the model, in which the longitudinal electric field 
does not exist near t he neutron-s tar surface (the model with 
free particle escape. [Arondl 19791 ; .Mestel 1999) . 

The evolutionary track of a neutron star in the P — sin x 
diagram does not depend on the magnetic field strength di- 
rectly. At that, the rate at which a neutron star moves along 
its evolutionary track depends on the magnetic field. Be- 
sides, the magnetic field strength affects the death line and 
the line of the transition to the propeller stage, so it can 
affect the evolutionary track through the track lengths of 
a neutron star at different stages of evolution. Further, we 
want to note that the boundary of the region of orthogonal 
rotators does not depend on the magnetic field either. At 
last, we should emphasize that the bias of the track in the 
region of active radio pulsars is defined solely by the ratio 
((9| for the current losses model. Therefore, while construct- 
ing the evolutionary tracks by the initial parameters (direct 
way), one should keep in mind that Eq. ^ is very suscep- 
tible to the small changes in the initial angle value in the 
case of the short (millisecond) initial period. 



Many of the newly discovered types of isolated neutron 
stars have long spin periods and may belong to the next evo- 
lutionary stages such as extinct radio pulsars or supersonic 
propellers, or even to orthogonal or aligned rotators. In that 
case, in order to estimate the age of the star we should re- 
construct the whole evolutionary track of a neutron star and 
consider a neutron star age as a sum of time intervals during 
which a star lives at each of the stages. 

As example, in this section we present reconstructions 
of evolutionary tracks for two NSs belonging to different 
classes: the Magnificent seven and anomalous X-ray pulsars. 
The choice is motivated at first by the fact that for both of 
them some information about inclinations angles is avail- 
able. For one of these NSs (RX J0720.4-3125) the angle is 
small or intermediate, as it appears from modeling it's X-ray 
light curve (see below). For another, AXP 1E2259-I-586, the 
angle is large, as it is clear f rom the double peaked X-r ay 
pulse profile (see, for example. I Woods fc Thompson! l|2004f )). 
Then both of these NSs has relatively high magnetic fields, 
which means that the rate of their magneto-rotational evo- 
lution is high, too. Finally, poor, but independent, age es- 
timates are available for them. In particular, the charac- 
teristic agefor RX J0720.4-3125 is r = P/2P = 1.9 Myr 
iHaberfcOOd l is larger than it can be esti mated from coolin g 
curves. Cooling age should be < 1 Myr (iPage et al.ll2004l ). 
In each of the following plots an independent age estimate is 
marked by an asterisk. The independent age estimates are 
taken as 1 Myr for RX J0720.4-3125 (excluding two plots 
with the characteristic ages less than the independent ones, 
in which cases the mark is placed at the origin, indicating 
that initial period and angle should be small) and 0.01 Myr 
for IE 2259-f 586. The latter is based on the as sociations of 
the s ource with supernova remnant CTB 109 l|Kaspi et al.l 
l2003h . 

At first, we consider the evolution of the X-ray dim 
isolated neutron star RX J0720. 4-3125 within the model of 
current losses, and reconstruct possible evolutionary tracks 
in P - sinx plane assuming that the angle between spin 
and magnetic axes at t he present time is Xobs ~ 50° 
(Perez- A zorin et al ] |2006l ). Basing on the present day period 
Fobs ~ 8.39 s and the angle Xobs, this object can be either 
in the region of active radio pulsars or in the region of ex- 
tinct radio pulsars depending on the magnetic field strength. 
We consider the reverse track reconstruction, i.e. knowing 
the parameters (the spin period and the inclination angle) 
of the neutron star at the present time we reconstruct the 
evolutionary track. After that we can find, how the initial 
parameters are related to each other. The initial spin period 
Po = 20 ms is used only for the age estimate of the neutron 
star. 

Here we reconstruct the evolutionary tracks of 
RX J0720. 4-3125 assuming the following three values of the 
magnetic field: B = lO" G, B = lO"'^ G, and B = 10^^ G. 
Knowing the magnetic field strength, we can define the 
death line of radio pulsars and the lines of the transitions to 
the stages of orthogonal rotators and supersonic propellers. 
After that we can determine if the point with the presently 
observed spin period Pohs and the presently observed incli- 
nation angle Xobs is located in the region of active radio 
pulsars or it falls to the extinct radio pulsar stage. Accord- 
ing to the position of this point on P - sin x diagram, we can 
define the next steps to reconstruct the evolutionary path. 
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Figure 1. The reconstructed evolutionary track of XDIN 
J0720.4-3125. B = 10^^ G, Fobs = 8.39 s, Xobs = 50°. Model 
based on current losses is used. The figure is valid for both hin- 
dered and free particle escape. The presently observed spin pe- 
riod and inclination angle correspond to the region of active radio 
pulsars. If the initial spin period Pq = 20 ms, the age estimate is 
r ~ 0.4 Myr. The "star" marks the initial parameters correspon- 
dent to the independent age estimate Ti„^=0.4 Myr. 



If we assume the magnetic field B — 10^* G, the present 
parameters of J0720. 4-3125 correspond to the point in the 
region of active radio pulsars. Recall, that during this stage 
the ratio sinx/P is invariant. Therefore, the evolutionary 
track of RX J0720. 4-3125 can be represented by a straight 
line between the points (Pobs; sin Xobs) and (0;0) (see Fig.[T|. 
The point defining the initial spin period Po and the initial 
inclination angle xo is located on this line; thus, if we know 
the initial spin period, we can unambiguously define the ini- 
tial angle between the magnetic and the spin axes. If the 
initial spin period is short enough (Po = 20 ms), the age of 
RX J0720.4-3125 can be estimated as 0.41 Myr HIJ- 

For lower magnetic field strengths the observed param- 
eters of RX J0720. 4-3125 correspond to the region of extinct 
radio pulsars. In that case, the evolutionary tracks of the iso- 
lated neutron star will be different for the following two mod- 
els of an accelerat ion region: the model wit h hindered par- 
ticle escape ( R uderman fc Sutherland 19751) and t he model 
with free particle escape (|Aronsll 19791 : lMestel|[l999l ) from the 
neutron star surface. 

Within the model with hindered particle escape plasma 
fills only the inner region of the pulsar magnetosphere, and 
a neutron star slows down according to the magnetodipole 
formula ([2)), therefore cosx/P is invariant for extinct radio 
pulsars. Hence, the track for this XDIN in the region of 
extinct radio pulsars in P — sin x plane can be represented 
as follows: 



'X = (1 - -P^ (cOSXobs/Pobs)^)° 



(18) 



In order to continue the evolutionary track of the star 
to the region of active radio pulsars we should find the pa- 
rameters, at which neutron star passes through the death 
line of radio pulsars (HSJ. As a result, for RX J0720.4-3125 
we have: 



r>0.75 

i>i2 



/ cos Xobs \ ' 
V Pobs / 



0.35_B?2'''' s. 



(19) 



dO.75 
cos Xdl = Pl2 



f H^i^V '' ^ 0.027B? 

V Pobs / 



(20) 



Knowing the coordinates of the point, where the evo- 
lutionary track for the extinct radio pulsar region intersects 
the death line, we can continue the evolutionary track of 
RX 0720.04-3125 in the region of active radio pulsars (the 
procedure is the same as it is shown for the example with 
B = 10" G). 

According to by the time of observation the neutron 
star lives as an extinct radio pulsar: 



T2 = 



dP 



2.4 X lO-ies^^p-i (1 - cos2 Xobs(P2)/(P„\J) 



s.(21) 



Finally, for the lifetime of a neutron star as an extinct 
radio pulsar at the time of observation we have: 



r2 = 2 X 10^^Bi2^-Po^bs cos ^ Xobs X 



(22) 



X In 



P? - -Pobs/cOS^ Xobs 



p2 

-^obs 



■Pobs/C0s2 Xobs 



Thus, the age of an extinct radio pulsar can be esti- 
mated as: 

r = n + T2, (23) 

where ri is the time neutron star lives as radio pulsar (|ll|l 
and T2 is a lifetime in the region of extinct radio pulsars. 

Fig. [2] shows the reconstructed evolutionary track for 
RX J0720.4-3125 within the model with hindered parti- 
cle escape from the neutron star surface for magnetic field 
strengths B = lO" '^ G. 

If the model with free particle escape from the neutron 
star surface is valid then an angle between the spin and the 
magnetic axes tends to 90° for both active and extinct radio 
pulsar regions. So, we can simply follow the same procedure 
for the extinct radio pulsars as we define for the active radio 
pulsars (see Fig.[T]and the example with B = 10^* G above) 
and use (|lip . Assuming the spin period of the star at birth 
was Po = 20 ms, the age of X-ray dim isolated neutron star 
RX J0720. 4-3125 within the model with free particle escape 
can be estimated as 11 Myr and 2 Myr for B — 10^^ G and 
B — lO^'^'^ G, respectively. It should be noted here, that the 
theory of the energy release associated with the longitudi- 
nal currents in t he neutron star mag netosphere was gener- 
ally developed bv lBeskin et al.l (|l993l) under the assumption 
of th e hindered particle escape ( Ruderman fc Sutherland! 
Il975h . Thus, if particles leave the surface of a neutron star 
freely (| Arons| [T979 : Mcstol 1999), expression i^U^ is a rough 
approximation of the period derivative and so the ages es- 
timated under the assumption of the free particle escape 
can be altered by a factor of a few for the more accurate 
spin-down rate expression. 

In addition, in Figs. [3H4] we present the evolution- 
ary tracks of AXP 1E2259+586 for a set of possible val- 
ues of magnetic field strength and presently observed an- 
gle between the magnetic and the spin axes (see data 
on AXPs in genera l and 1E2259-I-586 in particular in 
IWoods fc ThompsonI (|2004l )). Finally, in Figs.[5H6]we recon- 
struct the e volutionary tracks o f RX 0720.04-3125, assuming 
Xobs = 5° (|Haberl at al.l 120061 ). The validity of the model 
with hindered particle escape from the neutron star surface 
is assumed in these figures. 
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Figure 2. The reconstructed evolutionary track of XDIN 
J0720.4-3125. B = lO^^-S g, Pobs = 8.39 s, Xobs = 50°. Model 
based on current losses is used. The figure is valid for hindered 
particle escape. The presently observed spin period and inclina- 
tion angle correspond to the region of extinct radio pulsars. If 
the initial spin period Pq = 20 ms, the lifetime at the active 
radio pulsar stage is ri fa 1.5 Myr, and the lifetime at the ex- 
tinct radio pulsar stage is t2 ~ 4.3 Myr. The total age estimate 
is r ft; 5.8 Myr. The asterisk marks the initial parameters corre- 
spondent to the independent age estimate Tin(j=l Myr. 
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Figure 3. The reconstructed evolutionary track of AXP 
1E2259-I-586. B = 6 x lO" G, Pobs = 6.98 s, Xobs = 80°. Model 
based on current losses is used. The figure is valid for hindered 

particle escape. The presently observed spin period and inclina- 
tion angle correspond to the region of extinct radio pulsars. If 
the initial spin period Po = 20 ms, the lifetime at the active 
radio pulsar stage is ri fa 0.8 Myr, and the lifetime at the ex- 
tinct radio pulsar stage is T2 ~ 0.5 Myr. The total age estimate 
is r Ri 1.3 Myr. The asterisk marks the initial parameters corre- 
spondent to the independent age estimate rind=0.01 Myr. 



As, most probably, both objects left the region of ac- 
tive radio pulsars, then the track reconstruction shows that 
initial periods of these objects are long enough. 
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Figure 4. The reconstructed evolutionary track of AXP 
1E2259+586. B = 6 x lO" G, Pobs = 6.98 s, Xobs = 85°. Model 
based on current losses is used. The figure is valid for hindered 
particle escape. The presently observed spin period and inclina- 
tion angle correspond to the region of extinct radio pulsars. If 
the initial spin period Pq = 20 ms, the lifetime at the active 
radio pulsar stage is ti Ri 0.7 Myr, and the lifetime at the ex- 
tinct radio pulsar stage is t2 ~ 0.6 Myr. The total age estimate 
is T 1.3 Myr. The asterisk marks the initial parameters corre- 
spondent to the independent age estimate Tind=0.01 Myr. 
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Figure 5. The reconstructed evolutionary track of XDIN 
J0720.4-3125. B = lO^^.s g, p^^^ = 8.39 s, Xobs = 5°. Model 
based on current losses is used. The figure is valid for hindered 
particle escape. The presently observed spin period and inclina- 
tion angle correspond to the region of extinct radio pulsars. If the 
initial spin period Pq = 20 ms, the lifetime at the radio pulsar 
stage is ri 1.4 Myr, and the lifetime at the extinct radio pulsar 
stage is T2 Ri 19.3 Myr. The total age estimate is r r; 20.7 Myr. 
The asterisk marks the initial parameters correspondent to the 
independent age estimate Tind=l Myr. 

4 DISCUSSION 

In this section we provide a brief discussion on several topics 
related to our analysis. 

Wc discuss possible modifications of models due to non- 
electromagnetic effects. Then we discuss the possibility to 
choose between free and hindered paxticle escape basing on 
the data on inclinations angles of extinct pulsars. 
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Figure 6. The reconstructed evolutionary track of XDIN 
J0720.4-3125. B = lO" G, Fobs = 8.39 s, Xobs = 5°. Model based 
on current losses is used. The figure is valid for both hindered and 
free particle escape. The presently observed spin period and incli- 
nation angle correspond to the region of active radio pulsars. If the 
initial spin period Pq = 20 ms, the age estimate is t 0.34 Myr. 
The asterisk marks the initial parameters correspondent to the 
independent age estimate Tin;j=0.34 Myr. 

Then, as the main impact of poor knowledge of x evo- 
lution is on age estimates and population synthesis models, 
we touch these subjects. 

4.1 The role of invariants and modifications 

One of highlights of this paper can be summarized as fol- 
lows. If a neutron star is born rapidly rotating, it reaches 
the stage of low energy losses with very short spin period, 
unless the initial inclination angle of a star is very close to 
the value correspondent to maximal spin down. This con- 
clusion crucially depends on the existence of the invariants 
((21 or . The existence of the invariants is based on several 
simplifying assumptions: 

• A neutron star is an ideal sphere; 

• Magnetifc field is due to an ideal dipole; 

• Magnetic dipole braking or current losses are the only 
electro-magnetic processes which slows down the star; 

• There are no additional non electro-magnetic processes 
related to internal structure of a NS or to external agents. 

Clearly, all four assumptions can be violated. In such a case 
most of the results presented above have to be reconsidered 
taking into account additional processes or/and conditions. 

In particul ar, we want to pay attention to the old paper 
bv lMacvl(|l974l '). This author notes, that if the dominant in- 
ternal field is poloidal, then due to non-sphericity of a star 
the alignment process is working. If, oppositely, the toroidal 
component of the field dominates - then the counteralign- 
ment takes place. This arguments are particularly interest- 
ing as magnetars, especially SGRs, are supposed to have 
huge toroidal fields, which are dissipated on th e spin-down 
time scale of these objects and deform the star (IStella et al.l 
l2005l ). 

Another poss ibility is related to remnant discs 
iMenou et al.|[200ll ) . They can provide additional spin-down 



without significant changes in the inclination angle, espe- 
cially at early stages of evolution. So, when the disc becomes 
inefficient a pulsar already has a relatively long period, but 
did not significantly change it's inclination angle, yet. 

Also, a process of vortex line migra tion during spin- 
down can be important (|Rudermanl l200d ) . Effectively, this 
process results in counteralignment, and so can partly com- 
pensate alignment in magneto-dipole model. However, as 
alignment goes as cos x oc P and vortex line migration leads 
to counteralignment as sinx oc P", < n < 0.5, compensa- 
tion is only partial. 

Finally, for a rigid and non-spherical neutron star the 
free nutation process plays significant role regulating the 
tendency of th e magnetic axis to align of counter-align with 
the spin axis. (|Goldreichlll"970h showed that in case of mag- 
netic dipole braking model the nutation amplitude decreases 
if the inclination angle of a neutron star is less than ~ 55° 
and increases if the inclination angle is larger than ~ 55°. 

Electromagnetic processes can also lead to small 
changes in the angle x- Probably, it can be th e case for PSR 
B1931-I-24 which shows a peculiar behavior (|Kramer et al.l 
20065, if the explanation is in switching between m agneto- 
dipole and current losses (jBeskin fc Nokhrinall20o3 ). 

4.2 Free vs. hindered particle escape 

Potentially, in the case of the longitudinal current losses the 
distribution in the inclination angle of the neutron stars be- 
hind the death line can provide necessary information to 
choose between free ( Arons 1979; Mcstcl 1999) and hindered 
(|Ruderman fc Sutherlandll975f ) particle escape. The reason- 
ing is simple: if the model of free escape is realized, than 
(especially among the strongly magnetized neutron stars) 
most of the neutron stars should have x close to 90° . In the 
opposite case, there should be no neutron stars very close 
to X = 90° , but this does not mean that most of them have 
X ~ 0°. However, this easy approach can be done only if neu- 
tron stars are significantly evolved, i.e. if their present-day 
spin periods are much larger than the initial ones. Unfortu- 
nately, our analysis presented above shows that for the radio 
pulsars and magnetars the initial periods cannot be much 
smaller than the observed ones if either the magnetodipole 
or current losses model is valid. 

Among all highly magnetized NSs the Magnificent seven 
can have relatively small initial periods in comparison with 
the present day values, as it comes out after comparison of 
there age estimates based on the standard assumption and 
its analogue for current losses with independent estimates. 
If this is the case, then it is interesting to note, that they do 
not show overabundance of orthogonal rotatotors. This can 
be considered as a (weak) argument in favor of the hindered 
escape model. 

4.3 Age estimates for highly magnetized NSs and 
their initial spin periods 

For the current losses model we compare the standard char- 
acteristic ages and the age estimates with the account of 
the angle evolution for the different types of isolated neu- 
tron stars under the specific assumption that all considered 
isolated neutron stars still remain in the region of active ra- 
dio pulsars, i.e. they did not pass either through the death 
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line or through the boundary of the orthogonal rotator stage. 
Table 3 presents the results for one individual neutron star 
from each of the following types: HB-PSRs, SGRs, AXPs, X- 
ray dim isolated neutron stars and RRATs. For some sources 
the validity of the above assumption, of course, can be ques- 
tioned. But for the considered objects this assumption re- 
sults in the lower age estimate, and that is what we want to 
obtain. 

As there is no exact information available about the 
angle between the spin and the magnetic axes for almost all 
of these isolated neutron stars, at first we consider the low 
estimates of their ages assuming constant inclination angle 
equal to maximal losses position, r^. 

To get the lower age estimate for the case when sources 
are already at the stage of extinct pulsars, r^max, in the 
framework of the current losses model, we allow the angle 
X to evolve and find the maximum angle between the mag- 
netic and the spin axes at which the neutron star with the 
presently observed spin period could still be located in the 
region of active radio pulsars in the plane P — sinx (see 
Fig- 7)0 The fact, that for extinct pulsars in the scope of free 
particle escape model such an estimate provides the lowest 
possible age is obviousQ For the observed period an extinct 
pulsar has x > Xmax, so for the same initial period the one 
which is now extinct has larger also initial inclination angle. 
Then this means, that for all the same evolutionary param- 
eters it always had larger angle, and therefore, smaller P. 
To reach the same observed period from the same initial it 
takes longer for the one which is now extinct. For the hin- 
dered particle escape model we made numerical estimates, 
which demonstrate that even in these case for all considered 
parameters, which have been chosen according to parame- 
ters of the sources in Table 3, our value r,^max still is a lower 
age estimate for extinct radio pulsars. 

The age of the NS with x = Xmax is calculated accord- 
ing to Eq. (11), where xo = Xmax -Po/-Ptr. Here Ptr - is the 
observed period. 

All objects presented in Table 3 belongs to classes of 
highly magnetized NSs with long spin periods. For them 
magneto-rotational evolution goes quicker, so any evolution 
of X on a spin-down time scale should be very pronounced. 
Even without precise information about inclination angles 
of highly magnetized NSs, observations allow us to say that 
they cover a wide range: from nearly orthogonal, to small x- 
From Tables 1 and 2 one can see, that NSs with B larger 
~ 10^^ born with short periods rapidly, in less than few tens 
of thousand years or even faster, become aligned or orthogo- 
nal rotators without significant spin-down, unless they were 
born with x extremely close to the critical value of maximal 
losses. As the latter does not seem very probable for all of 
these sources (for such long periods and strong fields the 
proximity to the critical values should be of order of min- 
utes), we can conclude that initial periods of highly magne- 
tized NSs should be close to the presently observed values. 

Note that this conclusion is in correspondence with re- 



^ Of course, for PSR J1734-3333 this estimate represent the upper 
age boundary. 

^ We do not consider the possibility that any of these NSs is 
on the stage of orthogonal rotators, as all of them demonstrate 
significant spin-down rates. 
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Figure 7. The maximum possible inclination angle Xmax can be 
defined as the angle at which the neutron star with the presently 
observed period crosses the death line or the line of the orthogonal 
rotators region. The death line I I12II and the boundary line of the 
region of orthogonal rotators I I15I I intersect at the point with spin 
period Pj^t « 0.35-6^2^® s. If the presently observed spin period 
is less than the period Pjnti the inclination angle Xmax is the 
angle at which the neutron star with the presently observed period 
crosses the line of the orthogonal rotators region. Otherwise, if 
the presently observed spin period is longer than P^t then the 
maximum inclination angle Xmax is determined by the death line. 

cent investigation of prope rties of supernova re mnants as- 
sociated with magnetars (|Vinc fc Kuipej l2006l l. The au- 
thors do not find evidence that these neutron stars were 
formed from rapidly rotating protoneutron stars. Thus, it 
is likely that magnetars, as active radio pulsars, can be 
born with comparat ively long periods. Recent calculation 
(jShibata et al.l |2006| ) support the conclusion that highly 
magnetized N Ss are born slowly rota ting. Furthermore, 
[Ferrario fc Wickramadinghd (|2003| . 120061 ) discussed that the 
neutron stars with strong magnetic fields tend to be born 
as slow rotators, so initial spin period of a neutron star may 
depend on the magnetic field strength. On the other hand, 
for magnetars non-electromagnetic contributions to the evo- 
lution of X can be important. 

4.4 Population synthesis 

One of areas in which poor knowledge of the evolution of x 
can have strong impact is population synthesis of binary and 
isolated NSs. Below we mainly focus on the case of solitary 
objects, but in many aspects this discussion can be applied 
to binaries, too. 

The inclination angle x explicitly appears in the for- 
mulas for spin-down rate and energy release in both models 
(magneto-dipole and current losses). The death line should 
also depend on this parameter, and in the case of the current 
losses model this is shown explicitly. Finally, extreme val- 
ues of this parameter can result in low energy losses. Taken 
all together, it appears that modeling the evolution of this 
parameter is at least as important as the evolution of the 
magnetic field. 

Most (if not all) of population synthesis models accept 
the standard assumption about spin-down, which is Eq.(l) 
with sinx = 1 =const. The death line for radio pulsars is 
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Table 3. Age estimates for different classes of isolated neutron stars in the scope of current losses model for the constant inclination 
angle X = 0° and for the inclination angle evolving to the maximum possible value Xmax • We assume that a neutron star is located in 
the region of active radio pulsars on the P ~ sin x diagram and the initial spin period is taken to be Pq = 20 ms. 



Name 


Class 


P, s 


B, G 


Tci, Myr 


Xmax, deg 


Txmax, Myr 


J1734-3333 


HB-PSR 


1.17 


1014 


0.05 


89.2 


0.1 


1806-20 


SGR 


7.5 


1014 


0.3 


77.5 


0.5 


1E2259+586 


AXP 


6.98 


6 X 10" 


0.6 


64.3 


0.8 


J0720.4-3125 


XDIN 


8.39 


1014 


0.3 


71.4 


0.5 


J1819-1458 


RRAT 


4.26 


5 X 10" 


0.5 


83.7 


0.9 



taken in the form of Eq. (4). It is often believed that the 
uncertainty in the initial values of x can be hidden in uncer- 
tainties of B. In some sense, this is true. But only when we 
deal with active radio pulsars, and if the inclination angle 
does not significantly evolve. On later evolutionary stages 
(propellers, accretors) x ^^id B "decouples". Such "decou- 
pling" means that we become interested in actual values of 
parameters, not in some particular combination. In particu- 
lar, to predict properties of isolated accreting NSs it is very 
important to know if most of them are expected to be aligned 
rotators (and so, non-pulsating sources), or, oppositely, they 
tend to be orthogonal rotators showing clear pulsations. 

Also, the death line should depend on both x a-nd B, 
and at least in the case of current losses the combination 
of these parameters in the definition of the line is different 
from the combination in the spin-down rate. The use of the 
standard assumption, which, as demonstrated above, does 
not have solid theoretical grounds, makes many results of 
population synthesis not very trustable. So, if the aim is 
to produce a kind of a global p opulation synthesis, like the 
one presented in (jPopov et al.l [2000) , including post-PSR 
stages, then it is necessary to follow the evolution of x and B 
separately. Unfortunately, it has never been done, and there 
are serious difficulties in fulfilling such a program, as none 
of theories gives such an evolution in good correspondence 
with observations of radio pulsars. 



5 SUMMARY 

We revisited magneto-rotational evolution of NSs with ac- 
count of the evolution of the angle x between spin and 
magnetic axes. In our consideration we used two models of 
energy losses: magneto-dipole and longitudinal currents. In 
general, the results obtained in the framework of the former 
are not new, but are presented for comparison and illus- 
tration. We compare output of our calculations with obser- 
vations, and with estimates according to the standard as- 
sumption widely used for age estimates and in population 
synthesis. In this assumption the spin-down is calculated 
due to magneto-dipole model with constant x = 90°. 

Our calculations demonstrate that for initial angles not 
too close to the critical ones (90° in the case of magneto- 
dipole losses, and 0° - for current losses) NSs do not be- 
come normal extinct pulsars crossing the death line, but 
appear in the region of low losses (aligned rotators in the 
case of magneto-dipole model, and orthogonal rotators in 
the case of current losses). At that, for small initial periods 
and strong (but still not magnetar) magnetic fields this evo- 
lution proceeds very rapidly, in tens or hundreds thousand 



years. This fact does not correspond with observations. We 
consider three possibilities to resolve this problem: 

• an essential part of radio pulsar are born with compar- 
atively long spin periods, close to the observed ones, or 

• the initial angles between the magnetic and the spin 
axes are close to critical values, or 

• both magnetodipole model and the model based on cur- 
rent losses require significant modifications and cannot be 
directly applied for detailed calculations of NS evolution. 

To test the first assumption it is necessary to perform 
a detailed population synthesis modeling. However, it seems 
possible, that the results of such computations would not 
support the hypothesis, as we expect an unobserved over- 
abundance of old pulsars with angles corresponding to low 
energy losses. Also, the total birth rate of pulsars should be- 
come higher, in contradiction with independent estimates. 

The second possibility that initial angles x are very close 
(within 1 degree) to one of the critical values seems too 
artificial, and does not have a solid physical ground. Clearly, 
it should be modified. 

Most probably, our understanding of NS evolution on 
the radio pulsar stage in not satisfactory, now, and so both 
theories cannot be directly applied for detailed computa- 
tions. Some non-electromagnetic effects can be important, 
especially during early stages of spin-down. They can be re- 
lated to the presence of debris discs around NSs, or with 
mechanical effects related to non-sphericity or to processes 
in NS interiors. These effects can be different for different 
types of NSs. 

Highly magnetized NSs seem to be most favorable 
sources to test predictions of theories of magneto- rotational 
evolution, and to probe initial parameters. The evolution of 
such stars is the most rapid, for them there are estimates of 
angles x. and independent ages estimates. Our calculations 
(including track reconstructions) show that if such sources 
as SGRs, AXPs, RRATs, and the Magnificent seven are be- 
yong the death line, then their initial spin periods should be 
close to the observed ones. This conclusion is valid if one of 
the discussed theories is working, or, in general, if the in- 
clination angles significantly evolve on the spin-down time 
scale. 

Poor understanding of the details of evolution related 
to changes of the inclination angle is mostly important for 
age estimates and for population synthesis studies. Taking in 
account important role of x in both models of energy losses 
(at least, equal to the role of the magnetic field strength) 
we can say, that neglect of its evolution is a "skeleton in 
a cupboard" of population synthesis models. The account 
of the evolution of x can be crucial for radio pulsar stud- 
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ies, and significantly impact calculations of the later stages 
(propellers, accretors, etc.). 

Summarizing, we can say that our understanding of NS 
magneto-rotational evolution is not enough to make detailed 
population models and to obtain exact ages. The use of 
the standard assumption of constant x contradicts present 
day models of radio pulsar energy losses and require serious 
modifications or basis. Most pr obably, de tailed numerical 
modeling ( Timokhin 2006; Co ntopoulos fc Spitkovsky . 20061 ) 
can provide the correct answer about details of magneto- 
rotational evolution, in the future, but, probably, additional 
inclusion of non-electromagnetic processes can be necessary. 
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